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Abstract 


Alternative food sources have become an important focus of research 
due to increased food demand coupled with reductions in traditional 
food productivity. In particular, substitutes for protein sources have been 
of increasing interest due to the unsustainability of traditional protein 
sources. Insects have been identified as a sustainable alternative to tradi- 
tional protein sources, as they are easy to produce and contain essential 
proteins, fats, and minerals. However, mass-rearing insects requires simi- 
lar considerations as farming traditional protein sources. To increase pro- 
ductively, growth and survival must be maximized at the highest possible 
densities while minimizing disease and food requirements. Here, we use 
the house cricket Acheta domesticus, a highly cultivated insect species, to 
investigate optimal densities for mass rearing at 14 days of age (4" instar). 
Nymphs were separated into density groups of 0.09, 0.19, 0.47, and 0.93 
cricket/cm? and monitored for growth and survival. Multiple regression re- 
vealed sex (p < 0.0001), density (p < 0.0001), and sex* density interaction 
(p = 0.0345) as predictors of growth rate. Survival to maturation was sig- 
nificantly reduced in both 0.47 (31%) and 0.93 (45%) cricket/cm* groups 
compared to the controls. A second experiment was then conducted to 
investigate the starvation resistance of adult crickets reared from 14 days 
of age at 0.09, 0.19, 0.93, and 1.86 cricket/cm?. A second multiple regres- 
sion analysis revealed only density (p < 0.0001) and to a lesser extent sex 
(p = 0.0005) to be predictors of starvation resistance. These results indicate 
that mass-rearing house crickets is most optimal at densities < 0.93 cricket/ 
cm?, where impacts on survival and starvation are minimal. Although these 
results have implications for cricket mass rearing, research on other end- 
points, including reproduction and the synergistic effects of other environ- 
mental factors, such as temperature and humidity, should be conducted. 
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Introduction 


Accelerating global population growth and consumption has 
put increasing strain on the world’s agricultural system (Godfray 
et al. 2010, Sorjonen et al. 2019). Our ability to produce enough 


food for the world’s growing population is also diminishing due 
to factors such as urban expansion, land degradation, climate 
change, and water scarcity (Sorjonen et al. 2019). As a result, lev- 
els of global undernourishment have risen to an all-time high of 
9.9%, with between 720 and 811 million people facing food inse- 
curity (United Nations 2020). Food demand worldwide is predict- 
ed to increase, with demand rising by 70% by 2070. As traditional 
food production methods are unlikely to fulfil current and future 
global food demands, alternative sources of food have been pos- 
ited (Sorjonen et al. 2019). 

One key area of food production that has received much at- 
tention due to its current unsustainability is traditional protein 
sources. Currently, global protein requirements are fulfilled largely 
by red meat, poultry, and seafood (Thavamani et al. 2020). These 
foods are a key source of several micronutrients, including iron, 
zinc, phosphorus, vitamin B6, and B12 (Ajwalia 2020). However, 
meat production is highly unsustainable, with major concerns 
surrounding environmental degradation, animal welfare, and 
the negative effects of excessive meat consumption (Thavamani 
et al. 2020). Along with plant-based foods, cultured meats, and 
mycoprotein-based foods, insect-based protein sources are grow- 
ing in popularity as meat/protein alternatives (Vandeweyer 2018, 
Thavamani et al. 2020). Insects, like traditional protein sources, 
provide high volumes of fat, protein, zinc, iron, and several key 
vitamins (DeFoliart 1991, Schabel 2010, Alexander et al. 2017, 
Thavamani et al. 2020). The rearing of insects for this purpose is 
also more efficient and sustainable than for traditional livestock, 
consuming orders of magnitude less water and producing signifi- 
cantly less greenhouse gas emissions (GHGs) (Lundy and Parrella 
2015). Insects are also superior sources of protein when compared 
to traditional protein products, containing a protein content of 
50-82% of the dry weight (Thavamani et al. 2020). In addition 
to producing significantly less GHGs, insects can utilize organic 
waste products (low-value diets) as food sources, which provides 
further prospects for sustainable insect rearing (Sorjonen et al. 
2020, Thavamani et al. 2020). Currently, the insect-as-food indus- 
try is expected to grow by 47% from 2019 to 2026, with 730,000 
tons being produced by 2030 (Savio et al. 2022). 
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Despite the increasing global production of insect food 
sources, there is a paucity of information and research on best 
farming practices, including optimal densities for farmed species 
(Hanboonsong et al. 2013). Optimizing rearing density is one of 
the main considerations when seeking to increase production while 
maintaining food safety. Currently, it is generally understood that 
the effects of overcrowding include alterations in insect behavior, 
physiology, and, most importantly, life history (growth and surviv- 
al) (Maciel-Vergara et al. 2021). However, industrial mass rearing 
densities are often unnaturally high, and the extent to which these 
densities affect insect behavior, performance, and well-being is not 
well understood (Francuski and Beukeboom 2020). 

In industrial insect rearing, mortality is a common negative ef- 
fect of increased density (Peters and Barbosa 1997). In the Mediter- 
ranean fruit fly, for example, densities of 1-15 eggs/cm? show in- 
creased mortality as density increases due to increased food compe- 
tition (Peters and Barbosa 1997). Insect overcrowding also increas- 
es mortality as a result of increased pathogen susceptibility caused 
by higher temperatures, reduced immune response, and nutrient 
deficiencies (Vergara et al. 2021, Savio et al. 2022). Viral, fungal, 
bacterial, and microsporidian pathogens are all frequently found to 
infect mass-reared insects (Vergara et al. 2021). Viruses such as Ache- 
ta domesticus densovirus (AdDV) can result in mass mortality and 
have even been behind the bankruptcy of a few cricket-rearing com- 
panies (Weissman et al. 2012). Increases in disease transmission 
ultimately result in the increased mortality of mass-reared insects. 

Insect growth rates are highly dependent on environmental 
conditions including temperature, predators, humidity, and 
resource availability (Barragan-Fonseca et al. 2018). Although 
most of these stimuli can be easily controlled, nutrient deficiencies 
due to often-overcrowded conditions in mass-rearing facilities 
can result in delayed maturation, survival, and reduced body size. 
These are highly negative outcomes when mass-rearing insects 
where high survival and increased body size are desired (Averill 
and Prokopy 1987, Agnew et al. 2002, Reiskind et al. 2004, Rivers 
and Dahlem 2014). In Anopheles gambiae (Jannat and Roitberg 
2013), larval mortality increased by 36% in high-density groups 
due to waste products decreasing substrate quality. In both the fruit 
fly Bactrocera tryoni (Morimoto et al. 2019) and Western tarnished 
plant bugs [Lygus hesperus (Brent 2010)], significant reductions in 
adult emergence, body weight, fecundity, and energetic reserves 
were evident in crowded treatments. Declines in growth and 
prolongation of larval developmental periods have also been 
shown in several species, including Culex quinquefasciatus (Ikeshoji 
and Mullai 1970), Trogoderma glabrum (Beck 1973), Tipula oleracea 
(Laughlin 1960), Endrosis sarcitrella (Andersen 1956), and Ptinus 
tectus (Peters and Barbosa 1997). Increased population density 
can also create greater social contact and irritation, increasing 
competition, aggressive encounters between individuals, fighting, 
and habitat destruction (Southwick 1971, Weaver and Mcfarlane 
1990). Mating can also be affected by overcrowded conditions 
(Gavrilets 2000). Martin and Hosken (2003) observed that 
females of the dung fly Sepsis cynipsea became less interested in 
remating when the population density increased. 

Interestingly, however, some species have shown a shortening 
of developmental time in response to increased density (Cohen 
1968, Hodjat 1969). In addition, shifts in development time and 
size have often been shown in both high-density (overcrowding) 
and low-density (group effects) treatments (Peters and Barbosa 
1997). As it is often difficult to determine what rearing densities 
are indicative of overcrowding in any particular insect, as well as 
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between-species behavioral and physiological differences, research 
on optimal densities in a variety of species is necessary. Investigat- 
ing the effects of rearing densities not only on survival and growth 
but also on resistance to starvation due to these overcrowded con- 
ditions (i.e., resource competition) is also of interest. 

Although 1900 different species of insects are consumed glob- 
ally, crickets are currently one of the most highly cultivated insect 
species, being utilized for both human consumption and as food 
sources for other species (Lundy and Parrella 2015, Francuski and 
Beukeboom 2020, Savio et al. 2022). In this study, we investigated 
optimal rearing densities in the cricket Acheta domesticus, consider- 
ing the negative outcomes of overcrowding on growth, survival, 
and starvation resistance in groups at densities relevant to mass 
insect rearing. Densities of 0.09, 0.19, 0.47, 0.93, and 1.86 cricket/ 
cm? were reared starting at 14 days of age (4" instar) and main- 
tained until death. Resistance to adult starvation was also investi- 
gated at these various densities to understand the impacts of vari- 
ous rearing densities on food competition, as this is likely to occur 
in high-density mass-rearing facilities. 


Methods 


Breeding colony.—House crickets A. domesticus were housed in 
a 93 x 64.2 x 46.6 cm acrylic terrarium covered with 1.5-cm thick 
Durofoam insulation. The colony crickets were from a homog- 
enous stock inbred for > 60 generations. Fans on top of the en- 
closed structure provided air circulation. Crickets were sustained 
at a 12 h-day—12 h-night photoperiod and at a constant tempera- 
ture of 29°C + 2°C using 60-volt UV heat lamps. Ad libitum food 
(Country Range MultiFowl Grower, 17% protein, Quick Feeds 
Feed Mill, Copetown, Canada) was provided. Ad libitum distilled 
water was made available in soaked cellulose sponges, and egg 
cartons were provided for shelter. Oviposition medium (Vigoro 
Organic Garden Soil, The Mosaic Co., Lake Forest, IL, U.S.A.) was 
present in small plastic containers (7 x 7 x 7 cm). Oviposition 
containers were collected after a 24-hour period and incubated at 
29°C + 2°C until eggs hatched (~ 14 days). 


Density variation.—A. domesticus were taken two weeks (14 days, 4" 
instar) post-hatch from a single colony oviposition container and 
randomly separated into four experimental groups. Crickets at this 
age/molt are approximately 0.25 mm in length. The experimental 
groups consisted of 50, 100, 250, and 500 crickets, creating densi- 
ties of approximately 0.09, 0.19, 0.47, and 0.93 cricket/cm?. re- 
spectively. Densities given are approximation. A range of densities 
below 0.93 cricket/cm2 was chosen to allow for density-dependent 
observations, as densities above this had been observed in the lab 
to result in mass die-offs. All groups were given a 2 x 2 egg carton, 
which lowered the density slightly. The 0.09 cricket/cm? group was 
chosen as the control, as the lowest-density group is expected to 
maximize resources per cricket and thus have fewer negative inter- 
actions. Experimental housing containers were 29 x 18.5 x 12 cm, 
and the crickets were housed for life and in the same conditions 
as the colony. Cricket containers were continuously monitored for 
cleanliness, as extremely dirty rearing environments may nega- 
tively impact growth and development. Containers were cleaned 
(crickets were moved into fresh containers) at a minimum of once 
a week; however, in the higher-density groups, cleaning occurred 
approximately every other day due to increased mortality and ex- 
crement production. Food and water were replaced daily to ensure 
the same resource availability in each group. 
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Life-history measurements.—Maturation in A. domesticus is denoted 
by the adult molt in which wings develop and individuals become 
sexually mature. Females are easily identifiable by the fully-devel- 
oped ovipositor. Crickets reach adulthood approximately 5-60 
days post-hatch. Sex, maturation mass, and development time 
(number of days post-hatch) were recorded for all density groups. 
Sample sizes for growth rate were 0.09 cricket/cm? (n = 41), 0.19 
cricket/cm? (n = 68), 0.47 cricket/cm? (n = 142), and 0.93 cricket/ 
cm? (n = 224). Sample sizes for the proportion of individuals sur- 
viving maturation were 0.09 cricket/cm? (n = 50), 0.19 cricket/ 
cm? (n = 100), 0.47 cricket/cm? (n = 250), and 0.93 cricket/cm? 
(n = 500). Mass was measured using an Accuris analytical balance 
with a readability of 0.001 g +0.002 g. Maturation mass (g) and 
development time were employed to calculate the growth rate. The 
number of crickets that matured was used to determine the pro- 
portion of individuals that successfully matured. 


Starvation resistance.—After determining the impacts of various 
rearing densities on growth and survival parameters, a second 
experiment was conducted to determine the potential responses 
to adult starvation at various rearing densities. Individuals from 
a second oviposition container were arranged as described above. 
Group densities were slightly altered based on our initial results 
that had suggested that although the proportion matured was 
reduced at densities > 0.47 cricket/cm?, growth rate was not af- 
fected. The experimental groups were separated and maintained 
as described above. Experimental groups included 50, 100, 500, 
and 1000 individuals, representing densities of approximately 
0.09, 0.19, 0.93, and 1.86 cricket/cm?, respectively. Three to four 
weeks post-maturation, 10 females and 10 males from each den- 
sity group were weighed and placed in individual containers to 
prevent cannibalism. The mass of each individual was recorded to 
determine whether increased body mass is related to starvation re- 
sistance. The number of surviving individuals in the 1.86 cricket/ 
cm? group consisted of 7 males and 10 females. For the starvation 
treatment, each individual cricket from each group was placed 
into a small cylindrical container and covered with plastic wrap 
secured by a rubber band. Holes were added to the plastic wrap to 
provide ventilation, and the sex and group of each individual were 
noted on the container. Ad libitum water was provided by plac- 
ing a water-soaked cellulose sponge in the container; the sponges 
were re-soaked daily. Although the crickets consumed the cellu- 
lose sponges, they provided insufficient nutrients. Mortality was 
noted daily and used to determine longevity. 


Statistics. —To determine the effect of sex, density, and possible in- 
teraction (sex*density) on both growth and starvation resistance, 
a multiple linear regression was conducted using the most appro- 
priate model. The best-fit model was determined using a step-wise 
AlCc comparison. A D’Agostino—Pearson omnibus normality test 
was conducted to ensure data was normally distributed. Survival 
curves were analyzed using a Gehan-Breslow-Wilcoxon survival 
analysis to determine differences in survivorship among density 
groups. To analyze differences in the proportion that survived to 
maturation, chi-square tests were applied. To determine significant 
differences between the various rearing densities and the control, 
a Fisher's exact test was applied to each rearing density compared 
to the control. Finally, significant differences in the mass of starva- 
tion groups were determined using a one-way ANOVA followed 
by Tukey’s multiple comparisons test. All statistical analyses were 
carried out using Prism Graph Pad 9. 
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Results 


Survival to maturation.—Chi-square tests indicated significant 
differences in proportion matured among the different rearing 
density groups (p < 0.0001). Fisher’s exact tests indicated 
significant differences between the 0.47 cricket/cm?(p = 0.0008) 
and 0.93 cricket/cm? (p < 0.0001) groups compared to the 0.09 
cricket/cm* density group. This constituted a 31% and 45% 
decrease in the proportion that matured, respectively (Fig. 1). 


Growth rate.—Growth rates were collected for all rearing density 
groups and are reported as mean growth rate + SD (Fig. 2). 
Prior to analysis, a D’Agostino-Pearson omnibus normality 
test was performed and confirmed normal distribution. We 
used AICc model selection to determine the best model for 
describing the relationship between sex, density, and growth. 
The best-fit model carried 77.46% of the cumulative model 
weight and included two predictors (sex and density) with 
interaction effects F (3, 471) = 70.79, p < 0.0001, R*= 0.3108; y 
= 8.277 + 1.18681 - 1.09482 — 0.576363. Results suggest that 
sex (p < 0.0001), density (p < 0.0001), and, to a lesser extent, 
sex* density (p = 0.0345) are significant predictors of growth 
rate. The multiple regression results are outlined in Table 2. 
A summary of growth and development measurements is 
outlined in Table 1. 


Mass of starvation groups.—The mass (g) of each male and fe- 
male A. domesticus used for starvation-resistant treatment was 
recorded immediately prior to experimentation and are rep- 
resented as mean mass + SD (Fig. 3). A one-way ANOVA indi- 
cated significant differences between groups F (7, 69) = 58.47, 
p < 0.0001, with a Tukey’s multiple comparison test indicating 
significantly reduced masses (p < 0.0001) in both the 0.93 and 
1.87 cricket/cm? females compared to the 0.09 cricket/cm? fe- 
male controls. Significant reductions in mass were also detected 
in 0.93 (p = 0.0018) and 1.87 (p < 0.0001) cricket/cm’* males 
compared to the 0.09 cricket/cm? male controls. Between-sex dif- 
ferences (p < 0.0001) were also detected in the 0.09, 0.19, and 
0.93 cricket/cm? groups. 


0.09 0.19 0.47 


cricket/cm? 


REKK 


0.93 


Fig. 1. The proportion of Acheta domesticus reaching maturation 
in each rearing density (0.09, 0.19, 0.47, and 0.93 cricket/cm? ). A 
chi-square test indicated significant differences in survival among 
rearing densities (p < 0.0001). A Fisher's exact test showed signifi- 
cant differences in proportion matured between 0.19 cricket/cm? 
(p = 0.0008) and 0.93 cricket/cm? (p < 0.0001) compared to the 
0.09 cricket/cm? density group. 


mM Matured 
G Non-Matured 
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Table 1. Summary of growth parameters of each Acheta domesticus rearing density group. 


Experimental Group N Growth rate Percentage of 

(mg/day) controls 

0.09 Male (control) 14 7.629 

0.09 Female (control) 27 8.941 

0.19 Male 34 8.239 8.00% 

0.19 Female 34 9.109 1.88% 

0.47 Male 71 7.830 2.63% 

0.47 Female 71 8.963 0.25% 

0.93 Male 112 7.234 - 5.18% 

0.93 Female 112 7.829 - 12.45% 

14 


Growth Rate (mg/day) 


Density (cricket/cm”) 


Fig. 2. Density-dependent effects of various rearing densities on 
the growth rate of male and female Acheta domesticus. Values are 
represented as mean growth rate of each group + SD. Growth rates 
were calculated by dividing the mass at maturation (mg) by the 
time taken to reach maturation (days) of each individual. Multi- 
ple regression analysis determined that sex (p < 0.0001), density 
(p < 0.0001), and to a lesser extent sex* density (p = 0.0345) inter- 
action were strong predictors of growth rate. 


Density—starvation.—Starvation resistance was measured as days 
survived after the removal of sufficient food for all density groups 
and is reported as mean survival + SD (Fig. 4). Prior to analysis, 
a D’Agostino-Pearson omnibus normality test confirmed normal 
distribution. We used AlCc model selection to determine the best 
model to describe the relationship between sex, density, and star- 
vation resistance. The best-fit model carried 72.65% of the cumu- 
lative model weight and included two predictor values: sex and 
density (F (2, 74) = 28.05, p < 0.0001, R? = 0.4312; y = 10.46- 
3.661 B1-2.70182). The results suggest that both sex (p = 0.0005) 
and density (p < 0.0001) are significant predictors of starvation 
resistance. Results of the multiple regression are summarized in 
Table 3. The Gehan—Breslow-Wilcoxon test showed significant dif- 
ferences in survivorship between the 0.93 cricket/cm? (p = 0.0030) 
and 1.86 cricket/cm? (p = 0.0008) groups compared to the lowest 
density (0.09 cricket/cm*) female group. Variation in survivorship 
was also evident in the 0.19 (p = 0.0328), 0.93 (p = 0.0063), and 
1.86 (p = 0.0001) cricket/cm? groups compared to the males in the 
lowest-density groups. 


Upper 95% CI Maximal growth Development Mass at 
rate time (Days) maturation (g) 
8.080 8.98 50279 0.388 
9.280 10.91 48.15 0.431 
8.518 10.00 50.00 0.413 
9.475 12.04 48.97 0.446 
8.026 9.68 49.51 0.387 
9.219 11.24 47.32 0.424 
7.392 9.87 48.21 0.349 
8.017 10.96 46.71 0.366 


Table 2. Multiple linear regression analysis with AlCc compari- 
son was used to determine the most correct model for predicting 
growth rate based on sex and density group. AlICc comparison was 
utilized to select the best model. Our model includes sex, density, 
and sex* density interactions, which carried 77.46% of the cumu- 
lative model weight. Each predictor value had a significant correla- 
tion with growth rate: sex (p < 0.0001), density (p < 0.0001), and 
sex* density (p = 0.0345). 


Variable _—_ Coefficient (B) SE 95% CI P Value 
Intercept S277 0.1398 8.002 to 8.552 <0.0001 
Sex 1.186 0.1885 0.815 to 1.556 <0.0001 
Density - 1.094 0.1991 -1.485 to-0.703 <0.0001 
Sex* Density - 0.576 0.2718 -1.110to-0.042 0.0345 


Table 3. Multiple linear regression analysis with AICc comparison 
was used to determine the most correct model for predicting sur- 
vival based on sex and density group. Our model includes both 
sex and density, with no sex*density interactions, which carried 
72.65% of the cumulative model weight. Each predictor value had 
a significant correlation with growth rate: sex (p = 0.0005), density 
(p < 0.0001). 


Variable _—_ Coefficient (B) SE 95% CI P Value 
Intercept 10.46 0.6495 9.164to 11.750 <0.0001 
Sex 2.70 0.7451 1.217 to 4.186 0.0005 
Density - 3.66 0.5405 -4.738to- 2.584 <0.0001 
Discussion 


Insects have been proven to be a valuable alternative source of 
protein, fat, and essential vitamins and minerals (DeFoliart 1991, 
Schabel 2010, Zaelor and Kitthawee 2018). Utilizing insects as food 
may fill the gaps in our ability to meet current and future food de- 
mands. Crickets such as A. domesticus are one of the most cultivated 
insect species globally (Lundy and Parrella 2015, Zaelor and Kit- 
thawee 2018, Francuski and Beukeboom 2020). Acheta domesticus is 
considered an excellent candidate for mass-rearing endeavors due 
to its low food requirements and beneficial nutritional profile (Fer- 
nandez-Cassi et al. 2019). To mass rear insects for these purposes, 
it is vital to understand optimal densities for increasing individual 
size and survival (output) while reducing disease and stress. In this 
study, we investigated the impacts of various rearing densities on 
A. domesticus life-history features as well as their ability to survive 
prolonged starvation. We reared experimental groups for life-history 
analysis (growth rate and survival to maturation) at 0.09, 0.19, 0.47, 
and 0.93 cricket/cm*. We determined that the number of individu- 
als reaching maturation was significantly reduced in the 0.47 and 
0.93 cricket/cm? density groups (Fig. 1). This represented a decline 


JOURNAL OF ORTHOPTERA RESEARCH 2023, 32(1) 


S. MAHAVIDANAGE, T.M. FUCIARELLI, X. LI, C.D. ROLLO 


ae 


Mass (g) 


Density (cricket/cm?) 


Fig. 3. Mass of adult Acheta domesticus used for starvation treat- 
ment from each density group. Values represent the mean mass of 
each group + SD. The mass of each individual was recorded prior 
to starvation treatment. A one-way ANOVA indicated significant 
differences between groups (F (7, 69) = 58.47, p < 0.0001), with a 
Tukey's multiple comparison test indicating significantly reduced 
masses (p < 0.0001) in both the 0.93 and 1.87 cricket/cm? females 
compared to the 0.09 cricket/cm? female controls. Significant 
reductions in mass were also detected in the 0.93 (p = 0.0018) 
and 1.87 (p < 0.0001) cricket/cm? males compared to the 0.09 
cricket/cm* male controls. Although not represented on the graph, 
between-sex differences (p < 0.0001) were also detected between 
males and females in the 0.09, 0.19, and 0.93 cricket/cm? groups. 


of 31% and 45% compared to the lowest density 0.09 cricket/cm? 
group but still resulted in 142 and 224 individuals maturing, re- 
spectively. For growth rate, multiple regression analysis found that 
both sex (p < 0.0001) and density (p < 0.0001) have highly signifi- 
cant predictive power for growth rate. Interaction between sex and 
density (p = 0.0345) also had a slightly significant impact on growth 
rate (Fig. 2). A decline in growth rate was evident between the high- 
est density males and females compared to their lowest density con- 
specifics, constituting a 5.18% and 12.45% decline, respectively. 
Although our study did not indicate large declines in growth 
due to increased density, our results suggest that density is a strong 
predictor of growth rate. Prior studies have indicated that both 
increased mortality and decreased growth are expected due to 
overcrowding (Peters and Barbosa 1997, Zaelor and Kitthawee 
2018). These impacts on survival and growth are likely due to in- 
creasing population density, which has been shown to increase 
competition, physical injury, and stress (Parry et al. 2017). In an 
early study on the rearing densities of larval American cockroach 
Periplaneta Americana (Wharton et al. 1967), increasing density 
reduced both survival and growth. Other studies show similar 
trends, although the magnitude of survival and growth reduc- 
tions seem to vary among even closely related species. A study 
by Parry et al. (2017) found that both survival and growth were 
significantly affected by density; the relationship was not linear 
and was significantly different between the five blow fly species 
used. Studies on the English grain aphid Sitobion avenae (Xing et 
al. 2021) revealed decreases in both the growth and survival of 
early-instar nymphs with increased population density. The effects 
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Fig. 4. Starvation resistance of adult Acheta domesticus reared at 
various densities. All density and control groups were separated 
at 14 days of age (4 instar) and maintained until adulthood. 
Three to four weeks post-maturation individuals were separated 
and deprived of sufficient food. The number of days survived 
was recorded for all crickets and represented as mean + SD. 
Multiple regression analysis suggests both sex (p = 0.0005) and 
density (p < 0.0001) to be significant predictors of starvation re- 
sistance. The Gehan-Breslow-Wilcoxon test showed significant 
differences in survivorship between the 0.93 cricket/cm? (p = 
0.0030) and 1.86 cricket/cm? (p = 0.0008) groups compared to 
the lowest density (0.09 cricket/cm*) female group. Variation in 
survivorship was also evident in the 0.19 (p = 0.0328), 0.93 (p 
= 0.0063), and 1.86 (p = 0.0001) cricket/cm? groups compared 
to the lowest-density males. 


of density also seem to be age dependent, with some species be- 
ing more resistant to density in later life stages. For example, adult 
density does not significantly affect survival or reproductive traits 
in C. homnivorax (Berkebile et al. 2006). 

Varying results in the literature indicate that population den- 
sity effects are highly complex and species-specific, with some spe- 
cies being more resistant to the negative effects of overcrowding 
than others (Xing et al. 2021). Synergistic effects, such as tempera- 
ture and other environmental factors, should also be considered, 
as they may influence the effects of density, as was shown in a 
study using Sitobion avenae (Xing et al. 2021). This variability high- 
lights the need for species-specific data on density impacts. It is 
also likely that a more substantial reduction in growth may have 
appeared at higher densities than used here. In addition, as shown 
in Fig. 1, it is likely that the use of even higher densities may not 
only further impact growth but also survival. 

The second key aspect of this study was to investigate the 
impact of rearing density on starvation resistance. Intermittent 
food shortages are common in both the wild and in high-density 
mass rearing environments where competition for food is high 
(Zhang et al. 2019). Lack of food over extended periods of time 
can affect the growth, survival, and reproduction of individuals 
within the population (Zhang et al. 2019). Our chosen regression 
model suggested that sex (p = 0.0005) and to a larger extent den- 
sity (p < 0.0001)—but not the interaction between the two—are 
significant predictors of starvation resistance (Fig. 4). Analysis of 
survivorship curves indicated that most groups showed significant 
differences in survivorship compared to their same-sex lowest- 
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density group. The survival rate of the males in the 1.86 cricket/ 
cm2 density group was less than the males in the 0.09 cricket/ 
cm2 control group. It is important to note, however, that the mass 
of individuals used in the starvation resistance experiments was 
significantly lower in the 1.83 cricket/cm? and 0.93 cricket/cm? 
groups for both males and females (Fig. 3). The 1.83 cricket/cm? 
density group was not measured for growth rate, but it is likely 
that at this extremely high density the growth rate would be re- 
duced given the significant reduction in mass observed before star- 
vation. Our results are in line with to the manner in which insects 
are able to resist starvation in conditions in which migration is not 
feasible. Under starvation conditions, insects will undergo physi- 
ological modifications that alter their metabolism to help them 
cope (Zhang et al. 2019). They will first metabolize blood sugar 
(trehalose) and then lipids (triglycerides) to improve hunger re- 
sistance (Zhang et al. 2019). It is therefore likely that the reduc- 
tion in mass observed in the 1.83 cricket/cm* and 0.93 cricket/cm? 
disadvantaged these individuals in terms of their ability to break 
down sugar and fat stores as effectively as larger individuals. 

Although not always considered, sex often plays a key role in 
stress-related impacts on life-history features. Our results for both 
growth and starvation resistance showed significant contributions 
of sex on both variables (Tables 2, 3). Increased density has been 
shown to not only affect life-history traits but also alter interac- 
tions between different sexes (Rull et al. 2012, Parry et al. 2017). 
This may lead to alterations in fecundity and fertility (Rull et al. 
2012). Sex differences are expected, as females of this species are 
typically larger than males (Lyn et al. 2012). For starvation resist- 
ance, sex differences were evident in the 1.86 and 0.93 cricket/cm? 
groups, with females being more sensitive to starvation than males 
(Fig. 4). However, most studies have found female insects to be 
generally more resistant to food-related stress, as females are typi- 
cally larger and therefore have more nutrient stores. For example, 
Gaskin et al. (2002) found male Ceratitis capitata to be more nega- 
tively impacted by increases in density than females. The research- 
ers postulated that this was due to increased aggression and be- 
havioral costs to mate successfully in males (Gaskin et al. 2002). A 
study conducted on five species of blow flies reared at various den- 
sities indicated that females survived longer than males across all 
species used (Parry et al. 2017). Higher mortality in males versus 
females due to density variation has been recorded in Lucilia (Parry 
et al. 2017) sericata and Ceratitis capitata (Gaskin et al. 2002), Males 
of Cochliomyia homnivorax (Berkebile et al. 2006, Pitti et al. 2011) 
tend to show increased mortality under several rearing conditions, 
including high density, protein rich diets, and high temperature . 
While it is unclear why females in our study were more sensitive to 
starvation, trade-offs between cell maintenance and repair, result- 
ing in aging and death, and the energetic costs of egg production 
have been documented in female insects (De Loof 2011). The fe- 
males used in this study were sexually mature adults, potentially 
making them less able to mitigate the impacts of starvation. 

These results have profound implications for insect farming 
in which productivity is often deterred by the increased competi- 
tion and stress associated with high density (Zaelor and Kitthawee 
2018). Our results indicate that optimal densities for the mass rear- 
ing of A. domesticus are likely to be < 0.93 cricket/cm?, as this mini- 
mizes reductions in growth and maintains adequate resistance to 
starvation in adulthood. In addition, although survival to matura- 
tion is significantly reduced at this density, the number of individ- 
uals that do survive, in this case 224, is significantly greater than at 
densities with < 250 individuals. Thus, as the goal of mass rearing 
is to produce the largest number of individuals at the maximum 
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body size while reducing stress, we believe that < 0.93 cricket/cm? is 
optimal. These results should guide future mass-rearing endeavors 
to optimize production while reducing mortality and other nega- 
tive effects of overcrowding. It is also recommended that future re- 
search focus on a diversity of endpoints, such as how density influ- 
ences reproductive output, immune responses, and survivorship in 
this species to further inform optimal rearing densities. Synergistic 
effects between density and other environmental factors such as 
temperature, humidity, etc., should also be investigated. 
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